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We investigate the pseudospin-dependent density-energy relation (whose differential with respect 
to energy is density of states) of a monolayer graphene under intense terahertz laser field by exactly 
solving the time-dependent Schrodinger equation with help of Floquet's theorem. We find that 
psedospin polarization can be induced by a circular polarized terahertz laser field. The psedospin 
polarization in K and K' valleys can be exactly opposite sign when the electron densities in these two 
calleys are equivalent. Further more, we find that the psedospin polarization can be manipulated 
by the strength, frequency, and especially the polarization orientation of the field. 
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After the success of isolating the one carbon atom thick film 
graphene, ^'^ which is the thinnest material throughout scien- 
tific history, scientists are exploiting this new field^i^i^i^iii^ 
both theoretically and experimentally although it has been 
first studied more than 50 years ago.— Such a single layer 
graphene system contains two distinguishable carbon atoms 
A and B in each unit cell (A and B sublattice). The bond- 
ing TT and antibonding tt* bands touch with each other near 
the Fermi energy at two independent, degenerate and chi- 
ral symmetric Dirac points K and K' of the first Brillouin 
zone. Near these two points, the unusual dispersion which 
is approximately linear can form conically shaped valleys. It 
is intresting that the electron behavior in these valleys are 
govern by relativistic Dirac equation and the charge carriers 
can be regarded as massless Dirac Fermions^i^ with a effec- 
tive speed of light vf ~ l/400c. However, one of the most 
important potential applications of graphene is the graphene- 
based electronic devices. For example, the valley degree of 
freedom,«lS sublattice degree of freedom,— and layer index in 
graphene bilayer*^^ are suggested to be used in the same way 
as the electron spin is used in spintronics devices.— i^^ It have 
been proposed recently that intense terahertz (THz) field is 
one possible tool for contralinng and manipulating electron 
spinj^i^^^ In this paper, we propose a sublattice polariza- 
tion generator by using circular polarized terahertz laser. 

We describe the sublattice degree of freedom as a 
''pseudospin"^ The effective-mass Hamiltonian of an ideal in- 
finite graphene can be written as: 



//o(k) — hVF[cfxTkx + <Jyky\ , 



(1) 



where the Fermi velocity «f ~ 8 x lO^m/s, r = -|-(— ) repre- 
sent the states in K (K') valley. Here a is the Pauli matrix 
of pseudospin, o^ = ±1 denoting states on A (B) sublattice 
i. e. "up" ("down") state of pseudospin. We describe a cir- 



cular polarized THz laser field E(t) — — =(xcosnt -|- y sin Hi) 
with period To = ^, the vector and scalar potentials can be 
chosen as A(t) = — jS-[xs\nQ.t — ycosQt) and = under 
Coulomb gauge. With the gauge covariant derivative opera- 
tor P — —ihV — eA(i), the time dependent Hamiltonian can 
be written as: 



H(k,t) = hvF 



o'xT ( kxA — ;=7^ sin Q,t J +ay I ky ;= — cos Q,t 



V2hn 



V2hQ. 



(2) 

In order to understand the effect of circular polarized field 
clearer, we also need to know the effect of linear polarized field 
for comparison. We have ELinear(t) ~ EoxcosQt, A(t) = 
— ^xsinQt, the according Hamiltonian is: 



i/L 



r(k, t) = hVf 



O'xT ( kx + — TT- sin Qt ) + ayky 



(3) 



We solve the Schrodinger equation exactly by using the 
Floquet's theorem /^'^^ 



i^*(r,t)=H(k,t)*(r,t) 



(4) 



The solution is *I's,k(r, t) 

Ref. 193, •i's.kii) 
eigenvalue with s 






<l?s,k(i)- According to 
= <|i,,k(t)e-^«=''''', where g,(k) is the 
= ± represent two branchs of soul- 



a periodic function which satisfy 



H(k,t)0,,k(i). 
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Expand by Fourier series 



tion, (^s,k(t) 
(js(k)]<^s,k(t) 

</'s,k(i) = ^l^ao 't'",k^'"^" ^ 'P's.k IS a two Component spinor 
(''^s'k ' "^s k )"^- Then, the eigenvalues and eigenfunctions for 
circular polarized field and linear polarized field can be deter- 
mined by 



J 



[nn-gs(k)]. 



T / Ti + l 



[nfl — qs(k)]4>'s,k + ra'ax 



in + l 



"0?,k^] + VF{oxTkx + cryky)(j)"^i^ = , 
^ <^"k ] + VF{(^xTkx + cryky)(p",i^ = 

I 



(5) 
(6) 



Here a = Ox i: icy-u 






It is ob- viouse that such equations have periodic structure with the 



integer number of fl. If we replace q3(k) by gs(k) + luj, I 
is an arbitrary integer, one can find that the equations are 
unchanged, both gs(k) and gs(k) + leu could be the eigen- 
values. For convenience, we choose —Q/2 < gs(k) < /2. 
Moreover, since _ff(k,t) is Hermitian, the eigenvalues must 
satisfy X]^gs(k) = ^ J^ ° Tr7if(k, t)df,— therefore we have 
q_(k) — — q+(k) in our calculation. The density of states 
(DOS) is given by 



p(ti,t2)= /dk ^ <E..,k(ti)$I,, 

•/ ^ Li 



:(fe 



(7) 



field at r = 0, To/8, To/4 . We also plot the parabohc n-E 
curve n — J-f-r"^)^ without THz field. It is obvious that n-E 
curve can be remarkably modified by THz field, and it (as 
well as DOS) vary with time as the time-dependent electric 
field modulus. We pointed out that there is no pseudospin 
polarization in linear field case, the reason the symmetry be- 
tween sublattice A and B. If we interchange the index of A 
and B and let ky -^ —ky in Eq. ((6)1, the equation is un- 
changed. Since the summation over k in Eq. ((8)1, one can find 
Pa,a{T,uj) = pb,b{T,uj). 



it is a 2 X 2 matrix in pseudospin space. In order to translate 
into energy space, one can let T = ti + t2 and i = ti — i2.— 
After Fourier transformation with respect to t, one can obtain; 



Pil,i2iT,Uj) 



J{n-m)nT 



s — it n.m^ — oo 

x5[iu - (g^k) - (n + m)n/2)] , (8) 



where ^1(^2) = A 01 B, and i?5j,^2(s;n, 7n; k) = 
•j'*" k {4'TkV- Then the electron density for pseudospin state 

eis 

rEp(T) 

n5 = (l/27r)/ du;p^,^{T,Lu) . (9) 







,K 



nf + nf (n^' 



Once the total electron density n 
^A + iT-B ) in valley K (K') is given, one can determine the 
time-dependent Fermi energy Ep (T) {Ep (T)) according Eq. 
((9]). We define the pseudospin polarization in valley K (K') 



PX + PS (P 



K' 



pT +P§') 



In order to avoid 



treating insignificant singularity in numerical calculation, we 
interchange the integration dto and dk in Eqs. (|8]) (O, then 
we can obtain the relation between electron density and Fermi 
energy without presenting DOS. 




FIG. 1: (color online) Electron density versus Fermi energy 
in K (and K') valley under a linear polarized THz field with 
Eo = 1 kV/cm and Q. = 2-n: THz at T = (soUd curve), 
To/8 (dotted-dashed curve) and ro/4 (dashed curve). They 
are in contrast with the dashed parabola when the THz field 
is absent (dotted curve). 

First, in Fig. [l] we show the density-energy (n-E) curve 
whose differential with energy is DOS for linear polarized THz 
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FIG. 2: (color online) Electron density and pseudospin po- 
larization versus Fermi energy in K and K' valleys under a 
circular polarized THz field with Eq — 1 kV/cm and f7 = 27r 



THz. n'^ =n 
(dotted curve 



^ (solid curve) and P^ (dashed curve) — 
All the curves do not vary with time. 
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The circular polarized THz field case is much more inter- 
est. We plot both the n-E curve and pseudospin polarization- 
energy (P-E) curve of K and K' valleys in Fig. [2] Here, the 
DOS does not vary with time as the electric field modulus is 
a constant. One can see that the n-E curves of two valleys 
are exactly identical and the P-E curves are exactly oppsite in 
sign. This can be understood as the symmetry between t — + 
and — cases of Eq. ((Gjl . If we interchange the index of A and 
B for r = — case, and then let g3(k) — > 2nr2-|-gs(k), the eqau- 
tion will be the same as r = -I- case besides gs(k) —> — g3(k) = 
g_s(k). Because of the summation over s in Eq. (|8}, we can 

easily find that pA,Ai^) = Pb,b{'^) and pf,s(a)) = Pa,a{^)- 
However, the symmetry inside each valley we discussed of 
linear field case in above paragraph no longer come into ex- 
isting as the appearence of a^ and a^ in Eq. (|8}, therefore 

Pa,a(w) 7^ Pb,b{^) and Pa,a{'^) 7^ Pb.b{^)- A direct result is 
that once the given electron densities of two valleys that can 
be contraled by a gate voltage^^ are the same, in other words 
the Fermi energies of two valleys are the same, then the total 
pseudospin polarization P = P^ -\- P^ must be zero even 
though P'^ 7^ and P^ 7^ 0. We point out that one can use 
the so-called "valley filter " proposed by Recerz et al. in Ref. 
jlCJ ] to make the electron density (and Fermi energy) of one 
valley be very different from the other valley. Then, a total 
pseudospin polarization is obtained. 

Finally, we investigate the inffuence of the THz field pa- 
rameters such as field strength, frequency and polarization 
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FIG. 3: (color online) Electron density and pseudospin polar- 
ization versus Fermi energy in K valley with Eq = 0.5 (solid 
curve), 1 (dotted-dashed curve), 1.5 (dashed curve) and 2 
kV/cm (dotted curve), fi = 27r THz. 
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FIG. 4: (color online) Electron density and pseudospin polar- 
ization versus Fermi energy in K valley with E = 1 kV/cm and 
n = 1 (solid curve), 2 (dotted-dashed curve), 3 x 2-k (dashed 
curve)THz. An anti-orientated circular polarized fild case is 
also present (dotted curve). 



orientation in Fig. [3] and Fig. |4] We find that the n-E and 
P-E curves strongly depend on the field strength and fre- 
quency although. Especially, the polarization orientation can 
be utilized easily. The sign of pseudospin polarization can be 
changed if we apply an antiorientated THz field eventhough 
the n-E relation is unchanged. This feature may provide a 
very useful tool to manipulate the pseduspin polarization for 
graphene-based devices in the future. 



In Summary, we propose a method to generate and ma- 
nipulate pseudospin polarization in a monolayer graphene by 
using a circular polarized intense terahertz laser field. We 
solve the time-dependent Schrodinger equation with help of 
Floquet's theorem in this paper. It is pointed out that the 
pseudospin polarization in K and K' valleys can be exactly 
opposite sign when the electron densities in these two valleys 
are the same, and a total psedospin polarization can be ob- 
tained by utilizing valley filter which makes densities to be 
different in two valleys. Moreover, we find that the psedospin 
polarization can be manipulated by the strength, frequency, 
and especially the polarization orientation of the field. 
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